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Abstract 
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T) is evaluated within the leading ag In -r^- — approximation of QCD. In difference 
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from our earlier work on that subject, also the production of transversely polarized 
vector mesons is calculated. Special emphasis is placed on the role of the vector 
meson's qq light-cone wave function. In that context, conventional non-relativistic 
quarkonium models and a light-front QCD bound state calculation are critically 
examined and confronted with QCD expectations. Our numerical analysis finds a 
significant high momentum tail in the latter wave functions and a deviation from 
the expected asymptotic behavior of (pv{z,b = 0) oc z{l — z). We then design an 
interpolation to match the quarkonium models at large inter-quark separations with 
QCD expectations at small distances. We use these results to compare our predic- 
tions for the forward differential cross section of J/ijj photo- and electroproduction 
with recent experimental results from HERA. In addition, our earlier discussion 
of p° electroproduction is updated in light of recent experimental and theoretical 
enhancements. 

1 Introduction 

Diffractive vector meson production opens a precious window on the interface between 
perturbative QCD and hadronic physics. While elastic processes are commonly described 
through non-perturbative, phenomenological methods, as, for instance, soft Pomeron ex- 
change [m, hard inclusive reactions - most prominently deep inelastic lepton scattering 
- are, in a sense, exactly calculable as a consequence of the QCD factorization theorem. 
These two classes of processes now meet at HERA. However, similar to inclusive deep in- 
elastic scattering, also the amplitude for diffractive (coherent) production of vector mesons 
in deep inelastic lepton-nucleon scattering factorizes into a hard part calculable in pQCD 
convoluted with the non-perturbative off-diagonal gluon distribution in the target 0. A 
rigorous QCD-based proof of the factorization theorem for hard exclusive electroproduc- 
tion of vector mesons, valid to all orders in perturbation theory, was recently given in 
Ref. [^. This theorem holds if only short distances contribute, which is the case for the 
production of longitudinally polarized p° at sufficiently large Q^ or heavy flavor photo- 



and elect ropro duct ion |^. 

For large but non-asymptotic photon virtuality, the hard amplitude for exclusive vector 
meson production is sensitive to the transverse momentum distribution in the light-cone 
wave function of the qq leading Fock component of the produced vector meson [Q . This 
leads to a suppression of the asymptotic amplitude, i.e., to an interplay between the 
quark(antiquark) momentum distribution in the vector meson and the Q^ dependence of 
the corresponding cross section. That, in turn, allows to extract information on this wave 
function - and hence on the three dimensional distribution of color in the produced hadron 
- from the Q^ and the t dependences of the cross section. 

In this work, we focus the QCD analysis of Refs. and [^ on heavy quarkonium 
{J/ip and T) photo- and electroproduction. Furthermore, we extend the respective formal- 
ism, which in Refs. [0 and 0] was applied to the production of longitudinally polarized 
vector mesons only, to transverse polarizations as well. The important role the vector 
meson's qq light-cone wave function plays in diffractive photo- and electroproduction at 
non-asymptotic Q^ requires a detailed study of this quantity. Motivated by the large value 
of the quark mass in heavy quarkonia, we start from conventional non-relativistic potential 
models |^, ^, [7|, |^ and/or a non-relativistic light-front QCD bound state calculation [^. 
We then critically examine the respective wave functions and confront them with QCD 
expectations. 

In particular for the J /ip meson, our numerical analysis yields a significant value for the 
high momentum component in the respective non-relativistic wave functions (f)v{k)- For 
instance for the potential model of Ref. [^, the region - > 1, where the non-relativistic 
approximation is definitely inadequate, contributes over 30% to the integral Jd^k(j)y{k). 
Latter integral appears in the expression for the V -^ e^e~ decay width. This is illustrated 
in Fig. 6, and it is in line with the QCD prediction of large relativistic corrections to the 
corresponding bound state equations [|I^. Those large relativistic effects put the validity 
of a non-relativistic description of J /ip mesons - and, in particular, a non-relativistic 
evaluation of their production in high energy processes - seriously into question. Our 



analysis shows that the Q^ dependence of J /ip electropro duct ion, the photoproduction 
cross section ratios of T and J /ip mesons, and modifications in the t slope of those cross 
sections are good probes for the color distribution in the light-cone wave function of the 
vector mesons as well as the dependence of the parton distribution in the target on the 
produced meson's transverse size. In particular, these effects lead to an enhancement of 
the cross section ratio for diffractive electroproduction of T and J /if) mesons by a factor 
^ 10 for the same x as compared to the naive scaling estimate. This was discussed already 
inRef. §. 

In addition, if we express the non-relativistic wave functions in terms of light-cone 



coordinates, we find that they do not display the expected asjTiiptotic behavior |Tl 
J(fkt(pv L{z,kt) oc z{l — z) in the vicinity of 2 = or 2; = 1. This is illustrated in 
Fig. 8. Another mismatch between the non-relativistic and the light-cone approach ap- 
pears within the evaluation of the V -^ e+e^ decay width. When Vy^e+e- is calculated 
from the non-relativistic wave function (pyijx)^ a QCD correction factor, 1 — ^^^, appears 



1^ , which can be numerically large (^^^ ~ 0.5 for J/ip where we use as{J/ip) = 0.3) 



while no such term is present in the relation |T^ with the light-cone qq wave function 



(j)y{z,kt). This difference may be important in practice since the Schwinger formula for 
the positronium decay |T^ becomes inaccurate for charmonium where the high momentum 



component in the wave function is not small. To remedy these deficiencies, we designed 
an interpolation for the wave function of heavy quarkonia which smoothly matches the 
wave functions obtained at average inter-quark separations from non-relativistic potential 
models (or within a light-front QCD bound state calculation) with QCD predictions at 
small distances. 

The basic difference of the current work from Ref. [Q is that the formulae valid in 
leading order in „2i\^2 are derived by decomposing Feynman diagrams over the transverse 
distance between bare quarks, and that the quarkonium light-cone wave functions which 
we used respect QCD predictions for their high momentum tail. As for any hard process, 
the cross section is expressed through the distribution of bare quarks in the vector meson 



and not through the distribution of constituent quarks, as it has been assumed in Refs. 
p!5| and |l^. In the latter investigations, the cross section for diffractive photo- and 



electroproduction of J /ip mesons was evaluated in the BFKL approximation and while 
employing a non-relativistic constituent quark model. No corrections arising from the 



quark motion within the produced J/ip mesons were considered in Ref. [^. In a later 
work |]l6l , the authors then argued that the respective corrections are small within realistic 
charmonium models. This is at variance with our findings. In addition, our numerical 
analysis shows that the static approximation used in Refs. [IT3I and |TB[ is not in line 
with conventional charmonium models. Neglect of quark Fermi motion and related color 



screening effects in Refs. |T5| , |I6[ leads to factor ^ 3 suppression for the ratio of cross 
sections of photoproduction of T versus J/ip mesons as compared to the results in this 
paper. 

After outlining the basic formalism in Sect. 0, in Sect. ^, we discuss the heavy vector 
meson's light-cone wave function which describes its leading qq Fock state component. We 
then compare, in Sect. ^ with recent experimental results from HERA for J /ip photo- fl^ 



and electroproduction ||T8|. In Sect. ||, we update our discussion of p° electroproduction 
in light of recent experimental |T^ and theoretical [^ enhancements. We summarize and 
conclude in Sect. |[ 



2 The basic formalism 



2.1 The forward differential cross section 

In Ref. [0], the forward differential cross section for the production of longitudinally 
polarized vector mesons was deduced within the double logarithmic approximation, i.e., 
ttg In -7^ — In i ~ 1, with the result of 

QCD ^ 



da^^N^VN 



dt 



An^TvM^ 



t=o 



^chemQ 



ylVly 2 
^% 



a,{Q^){l+ip)xGN{x,Q^) 



Here, Ty stands for the decay width of the vector meson into an e~^e~ pair, P = j "^ ~ 
i 9inx ^^ ^^^ relative contribution of the amphtude's real part, and the leading twist 
correction 

2 Jdz Jd'^kt(pv{z,kt) 
accounts for the difference between the vector meson's decay into an e~^e~ pair and diffrac- 
tive vector meson production. Here (j)v{z,kt) is the wave function of the longitudinally 
polarized vector meson. We implicitly use the light-cone gauge which provides for an un- 
ambigous separation of the kf dependence of the meson (photon) wave function and gluon 
degrees of freedom. Note that the original formula deduced in Ref. lacks a factor of 4. 
This misprint has been corrected in Ref. [§]. In Ref. 0], it was shown that the formula 
in Eq. P) is valid also within the more conventional leading a.. In -r^ — approximation. 
Although, in principle, hard diffractive processes are expressed in terms of non-diagonal 
parton densities, an analysis [0] of the QCD evolution equations for non-diagonal parton 
distributions at small x shows that the difference between the diagonal and non-diagonal 
parton distributions is small in the kinematic region discussed in this paper. 

In Ref. [§], also next-to-leading order (NLO) as well as higher twist corrections were 
introduced. Firstly, it was argued that the strong coupling constant and the nucleon's 
gluon density have to be evaluated not at Q^ but at a Qlff- This is due to the so-called 
"rescaling of hard processes" which will be discussed in more depth later. And, secondly, 
a suppression factor T{Q'^) was deduced which measures the deviation of the cross section 
from its asymptotic prediction in Eq. (0), and which stems from the transverse Fermi 
motion of the quarks in the produced vector meson. This yields H 



(icr^£ 



N^VN 



dt 



A7^^TvM^ 



V^y-^V 2 rT.!r^2\ 



t=0 'ioiEAlQ 

with the correction factor 
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Tl'vT{Q') as{Qiff){l + i(3)xGN{x,Qiff) , (3) 



2 



T{Q' 
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r dz 
J z(l-z) 
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<Pviz, 


kt) 





(4) 



where 

M^^h) = --—^^ (5) 

^ + z{l-z) 

is the photon's qq hght-cone wave function, A^ is the transverse Laplacian A^ = J]{d/dki)'^, 
and where, for the production of hght mesons, the current quark mass was set to zero. 

In this investigation, we focus on the photo- and electroproduction of heavy quarko- 
nium {J/ip and T), and we extend the respective formahsm to the production of trans- 
versely polarized heavy vector mesons as well. Note that, for sufficiently heavy quark 
mass, the applicability of the QCD factorization theorem to the diffractive photoproduc- 
tion of transversely polarized vector mesons can be justified because the transverse size 
of quarkonium decreases with the mass of the heavy quarks. The result for the forward 
differential cross section for photo- and electroproduction of heavy vector mesons, which 
will be deduced in detail in the following, is 



do'^(*);^_,VN 



dt 



,.„ = ia^w^Lr 4T(Q^)hM„)(^+^0)^oA.,Qi„: 
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Here, riv is again the leading twist correction of Eq. (§), the factor T{Q'^)^ which was 
introduced in Ref. 0], accounts for effects related to the quark motion in the produced 
vector meson, and e = ^_'''~^2/2 is a parameter related to the (virtual) photon's polar- 
ization. Here, y is the energy fraction (in the target rest frame) transferred from the 
scattered lepton to the target. A value of e = corresponds to purely transverse polar- 
ization - which is always the case for real photons, i.e., Q^ = - and e = 1 refers to an 
equal mix of longitudinal and transverse polarizations. The latter is typical for HERA 
kinematics at large Q^ . The factor R{Q'^) parameterizes the relative contribution of the 
production of transversely polarized vector mesons as compared to the naive prediction, 
i.e., fj = R{Q^) ^ instead of simply ^ = ^. 

In difference from Ref.[0, the current quark mass was not set to zero, i.e., we kept 
leading powers over q2\^2 and not just ■^. This, in turn, yields for the correction factors 



r(g2) and i?(Q2): 



.2\2 



(Q^ + Am^y Jdz Jd^h Mz,h) A, (l)^iz,kt) 



m 



I 



dz 



z^il-z)' 



Jd'^kt (i)v{z,kt) At (py{z,kt) 



(7) 



4M2 Jdz Jd^kt (f)v{z, kt) At (j)^{z, kt) 

where we employed again 0^(-2, kt) of Eq. (|). 

The T{Q'^) and R{Q'^) displayed in the above constitute one of our main original 
new results. These formulae are derived by building a decomposition over the transverse 
distance between the bare quarks, i.e., over powers of qt:^^^- However, some caution 
is necessary at this point. The accuracy of this approximation for the calculation of 
R at Q^ ^ My can be questioned because of an enhancement of end point {z = and 
z = 1) contributions at large Q^. But in these kinematics, the production of longitudinally 
polarized vector meson would dominate [@. In order to be able to evaluate the correction 
factors of Eqs. (|^) and (|), we need the light-cone wave function of the qq leading Fock 
state in the vector meson. We will discuss this quantity in detail in the next section. 

Our master formula in Eq. (|^) yields a few fundamental predictions: 1) the cross 
sections raise with energy very rapidly due to the presence of the gluon density which 
increases fast at small x, 2) the t- slope is expected to be almost the same for all hard 
diffractive processes of the kind studied hereQ and 3) the production of longitudinally 
polarized vector mesons will dominate at large Q^. Note, also, that this is only a leading 
order analysis, and to achieve a non-ambiguous interpretation of the processes considered 
here it would be necessary to evaluate also more accurately NLO corrections as well as 
the higher twist the contribution of the \qqG) component in the light-cone wave functions 
of the photon and the produced vector meson. 



^This is because the i-slope in the hard vertex scales hke the maximum of the mass of the heavy quark 
and Q^. 



2.2 The color-dipole cross section 

As discussed at length in Refs. 0, 0] and [Q, due to the QCD factorization theorem and 
the large longitudinal coherence length, Ic ~ ^m x ^ associated with high energy (small 
x) diffractive processes, in leading order in a^ln-r^ — , the amplitude for hard diffractive 
vector meson production off a nucleon, depicted in Fig. 1, can be written as a product of 
three factors, 

A^i,)N^VN oc ^{-f* ^qq)-aggN-'^{qq-^V) (9) 

where \l/(7*^*^ -^ qq) is the light-cone wave function for a photon to split into a qq pair, 
(TggAr is the interaction cross section of the qq pair with the target nucleon, and \E'(gg -^ V) 
is the amplitude for the qq pair to transform into the vector meson V in the exit channel. 
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Figure 1: Feynman diagrams relevant for the evaluation of the amplitude for diffractive 
production of vector mesons, i.e., the 7*^*-' + A^ -^ V + A^ process, in leading a^ In -4^ — 
approximation. 



As was shown in Ref. @], for sufficiently large Q^ and longitudinal polarization, the 
above process is dominated by qq configurations where the quark and antiquark are sepa- 



rated by a small transverse distance b. Then, aqqN is the color-dipole cross section [^, |23 



(TggN{x,b) = ^b^ [asiQlff) xGN{x,Qlff)\ _q2+m2 2 _a ■ (^'^) 

Qualitatively, Eq. ([T0|) can be understood in the following way: The four diagrams of Fig. 
1 lead to an expression in the amplitude of the form 

(T-y*N oc I dz I (fit ... [2(l)^{z,kt) - (j)j{z,kt + lt) - (l)jiz,kt-lt)] , (11) 

where the Sudakov variable z denotes the fraction of the photon's momentum carried by 
one of the quarks, ±kt is their transverse momentum, and It is the gluons' transverse 
momentum. For small k, this yields 

a^*N oc dzcflt ... Il^t4>'y{z-,kt) . (12) 

Via Fourier transform into the transverse impact parameter space and after pulling out 
the wave function of the 7*, we obtain: 

aqqN oc 6^ . (13) 

The gluon density, xGn., arises as the diagrams in Fig. 1 represent not simple two-gluon 
exchange but rather the coupling to the full non-perturbative gluon ladder. For further 



details see Ref. ||2^ where the quantity aqqp^ was derived rigorously. In the following, we 
will show that - due to the large value of the current quark mass - the dominance of short 
distances holds for diffractive production of heavy flavors also for Q^ = and both for 
longitudinal as well as transverse polarizations. 

Note that, due to the difference in the invariant mass between the photon and the 
vector meson, the light-cone momentum fractions of the gluons in the initial and final 
state, (3i and /5/, are not the same, and therefore, in principle, an off-diagonal gluon 
distribution should enter into Eqs. (^ and (pn|). This was first recognized in Ref. [Q, 
and then elaborated on in Refs. pSl, p6i and ||21|. A simple kinematical consideration 



yields (3i ~ ^q'^+m'^ — -^ ^^"^ f^f ^ ^q^+m'^ ^ ^i where (/j) is the average transverse 
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moinentum of the exchanged gluons and Mj^ = {^ffzh) i^ ^^^ invariant squared mass 
of the produced qq pair. Within the as In j^ — approximation, the non-diagonal gluon 



QCD 



distribution is shown [^ to be not far - at the small x that are important experimentally 
- from the diagonal one. This is because, within this approximation, the appropriate 
energy denominators only weakly depend on /5j. 

2.3 Rescaling of hard processes 

As outlined in detail in Ref. 0, the parameter A, which fixes the scale in the gluon 
density and the strong coupling in Eqs. (H) and ([T0|), is determined by comparison with 
the longitudinal structure function, Fi(x, Q^) oc yG Niv, Q'^)\yr:^2 5x^ ^•^•' ^^ setting 

yGN{y,Q^) ^^ (x fd^b dz \(j)^*{z,b)\^ aggN{2.5x,b) , (14) 

where 

(l)^l{z,b) = 2Q zil - z) Ko (b^jQ-'zil - z) + m^^ (15) 

is the light-cone wave function of the qq leading Fock component in a longitudinally 
polarized virtual photon, and m is the current quark mass, which was set to zero when we 
evaluated Eq. ([T^) (since the contribution of charm quarks to Fl is small in the considered 
kinematics). In Eqs. [iVj and (|15]), b is the transverse distance between the quark and 



antiquark within the photon. The quantity A is adjusted such that the average b = b^^, 
which dominates the integral on the right hand side of Eq. (p!4D, is related to Q^ just 
via the equality 6^^ ~ ^- ^"^ other words, for the longitudinal structure function, the 
virtuality that corresponds to the dominant transverse distance b^^ is just the virtuality 
of the process. This yields A ~ 8.5 for x = 10^'^. 

In the same fashion, we can now rewrite the amplitude for diffractive vector meson 
production as 

A^*^N-.vN OC as{Qlff)xGN{x,Qlff) dzd% (f)^*^{z,b) b"^ (j)v{z,b) , (16) 

where we pulled the gluon density at an average b = by out of the integral, i.e., Qlff ~ ^, 
and with the qq leading Fock state light-cone wave function of the vector meson, 4>v{z, b). 
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In Fig. 2, we show by and Qlff for the fongitudinal structure function as well as for 
diffractive production of (longitudinally polarized) p°, J/ip and T mesons. The wave 
functions (j)v{z,b) that were used to evaluate Eq. (|16D will be discussed in more detail 
later. 
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Figure 2: Average transverse distances effective in the evaluation of the longitudinal 
structure function as well as for diffractive production of longitudinally polarized p°, J/ip 
and T mesons. Also shown are the resulting effective scales, Qlff, for diffractive vector 
meson production. 



It can be seen from Fig. 2 that the relevant transverse distances for p° electropro- 
duction are larger than those characteristic for the longitudinal structure function, i.e., 
bp{Q'^) > bcrj^iQ"^). Therefore, for p° production, the virtuality Qlff that enters in the 
argument of ^^(a;, Qlff) xGn{x, Qlff) is smaller than Q^. We find, to leading order. 



bviQ' 



which, for p° production, yields Qlff 



'^(TL\Qeff) 5 



(17) 



Q' 



XQ^) 



bv(Q^) 



. Our Eq. (^) is an approxi- 



mate relation designed to overcome the scale ambiguity which is inherent to leading order 
calculations. This "rescaling of hard processes" effectively relates the scales in different 
processes via the dominant qq distances in the respective quark loops. We termed this 
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"Q^ rescaling" in Ref. |Q. The difference between Qljf and Q^ indicates that substan- 
tial next-to-leading order corrections should be present in those processes. Applying the 
same method to J/ip and T production yields a Qlff which is significantly larger than 
the estimate Q = — ^-^ of Refs. [|1^, . Fig. 2 also indicates that the relevant trans- 
verse distances are small, and hence the QCD factorization theorem is applicable, for p° 
production at large Q^ and heavy meson photo- and electroproduction. 



2.4 Production of transversely polarized vector mesons 

The discussion in the above refers to the production of longitudinally polarized vector 
mesons only. For light vector mesons, the formalism at hand cannot be extended to trans- 
verse polarizations because of the endpoint singularities, i.e., the contribution from very 
asymmetric qq pairs with z ~ or 1, where non-perturbative effects dominate. For the 
production of heavy quarkonia QQ, when ^^(Mq) <^ 1 and -^^ ^ tt however, effects of 
large transverse distances are strongly suppressed. Here, tt is the radius of the hadron 
target. So, the production of transversely polarized heavy quarkonia can be legitimately 
evaluated using the QCD factorization theorem. At the same time, for Q^ ^ My the end 
point contribution (z ~ or 1) is enhanced in the amplitude for the diffractive elctro- 
production of transversely polarized vector mesons. Thus, the region of applicability of 
non-relativistic wave function models for heavy quarkonia to the production of transversely 
polarized vector mesons (but not of the applicability of the QCD factorization theorem) 
is restricted by the kinematical constraint Q^ < My. 

Employing the notations of Ref. [0], the wave functions of longitudinally and trans- 
versely polarized photons and heavy vector mesons can be expressed as 



/,AiA2 

7L 

/,AiA2 

7t 



2Qz{l-z)<p,iz,b)SX, (18) 

^^l^' = -2MvMz,b)6X , (20) 



13 



w' = 777^(^1 1 MzM^ll, (21) 



V. - 4137) I . I v^vi^,^;"A, , 



where 



0^(z,6) = iro(&\/QMl--2)+^') (22) 

and 0v(-2, &) refer to the qq hght-cone wave functions of the photon and the heavy vector 
meson, respectively. For the derivation of Eqs. (|20|) and (|2lD it was assumed, in hne with 
the non-relativistic character of heavy quarkonium, that, in the center of mass system, the 
vector meson's wave function is a pure angular momentum L = state. This selects spin 
S" = (or helicities A2 = — Ai) for the longitudinal polarization and S* = 1 (or A2 = Ai) 
for the transverse polarizations, with the same spatial wave function (f)v{z-,h). Here, Ai^2 
are the helicities of the quark and antiquark, respectively. For transverse polarization, 
the restriction through the wave function of heavy quarkonia selects the component in 
the wave function of the virtual photon which is proportional to the mass of the heavy 
quark. This is just opposite to the production of mesons built of light quarks where this 
component in the photon's wave function is negligible W% - 

This gives for the kernels of the longitudinal and transverse amplitudes: 

Vl{z, h) = \ Ea,a, <P':t' ^vl" = - ^Q^y ^(1 - ^) <^7(^. b) M^, b) , (23) 

2 

Vt(^, h) = i Ea,A. <P''r''^ ^Vr' = "^f^ U^^ &) M^, b) ■ (24) 



Note that in the limit 2 ~ | and My ^ 2m, Eqs. (^3]) and (|2^) yield the naive 
prediction f^ = (^) ~ ^^ for the production ratios of longitudinal to transverse po- 
larizations. Also, due to the non-relativistic ansatz for the vector meson's wave function, 
the spin structure of Eq. (^) is such that there is no azimuthal asymmetry. This is 



qualitatively different from diffractive two-jet production in deep inelastic scattering [^ 
Note, however, that in a fully relativistic description such an azimuthal asymmetry would 
appear also for diffractive production of transversely polarized vector mesons due to the 
admixture of a L = 2, S* = 1 component (to the standard L = 0, 5* = 1 state). Note, also, 
that the non-relativistic approximation to the light-cone wave function of transversely 
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polarized heavy quarkonia becomes questionable for the diffractive electropro duct ion in 
the limit Q^ ^ My. This is because, in this kinematics, the end point contributions 
2; = and z = 1 are enhanced. But in QCD, in variance from non-relativistic quarkonium 
models, the wave function at asymptotical Q^ should be such that Vr oc z{l — z). Such 
a behavior follows from the analysis of pQCD diagrams for the wave function of heavy 
quarkonia. This complication is practically unimportant because, in this kinematics, the 
production of longitudinally polarized heavy quarkonia dominates. 

Putting everything together, the factor T(Q^), which accounts for effects related to 
the quark motion in the produced vector meson, and the correction factor i?(Q^), which 
parameterizes the relative contribution of the transverse production, can be written in 
transverse impact parameter space as 

2 






(Q2 + Amy Jdz z(l - z) Jdb (Pviz, h) b^ (j)^{z, h) 



m 



4 /^0v(^,& = O) 

2 Jjrt^^!dbMz,b)b^<P,{z,b) 



4M^ Jdzz{l-z)Jdb(l)v{z,b)b^ (p^{z,b) 



(25) 
(26) 



where we used again (f)^{z,b) of Eq. (|22D. The T{Q^) and R{Q^), displayed in the above. 



are the leading expressions to order 



1 



02+4m2 



and they constitute our main original new 



results. They are related to the quantities given in Eqs. (^ and (H) simply via a two- 
dimensional Fourier transformation. 



2.5 Leading twist expressions and comparison with other kt- 
suppression estimates 

Note that the suppression factor T^Q"^) of Eq. ( p5[ ) and the transverse to longitudinal 
production ratio R{Q^) of Eq. (^^ have contributions from leading and non-leading 
twist. The corresponding leading twist expressions can be deduced by pulling the vector 
meson's wave function, (j)v{z,b), at 6 = out of the integral, i.e, by replacing (j)v{z,b) 
with (pviz, 0). The latter is equivalent to setting in the photon's wave function, (f)-y{z, kt) 
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of Eq. 



kt to zero after differentiation, and it yields 

/ \2 -.2 



TlAQ') 



J 2(1 



dz 



Q'^+Am? 



(1-^) U^+irf 



IT^iy 



0y(z,6 = O) 



/^M^,fe = o) 



(27) 



i?LT(Q') 



dz 



m 



2^(1-2)3 



^(1-z) 



0y(^,& = O) 



r_2 
J 2(1 



C(2 



(1-2) 



0v(z,6 = O) 



2(1-2) 



(28) 



This shows that, for these processes, a decomposition over twists is really an expansion in 
powers of 6^, and "leading twist" is equivalent to the 6 — > limit, i.e., to considering very 
small transverse distances (or "pointlike hadrons") only. Specific to heavy quarkonium 
production is that, in addition to neglecting kf/Q"^ and iv? /Q"^ terms as for light quarks, 
one also neglects terms of the form k'^/ra^. 

Note that the expressions (pTf ) and (p8|) are stringent QCD predictions for heavy quark 
production deduced in an expansion where m is considered as a large parameter. The 
leading term is proportional to the mass of the heavy quark in difference from light quark 
production where the leading term is proportional to the quark's transverse momentum. 
So, the formulae deduced in this paper cannot be smoothly interpolated to the limit of 
the zero quark mass. 

Furthermore, in the static limit oim ^ 00, which implies 0v'(^? ^t) = ^\Z — \) <Pv{kt) 
and My = 2m, the correction factors T{Q'^), R^Q"^), Tlt^Q"^) and Rlt{Q^) reduce to 

(29) 
(30) 
(31) 
(32) 



where 



T{Q^) 


-. 1 32 ^'^^ 

^ ^^Q2+4^2 


nm 


^ 1, 


TlAQ') 


^ 1, 


RlAQ') 


^ 1, 


(k^\ = 


_ Jd^ktk^Mkt) 



j'(Pkt (t)v{.kt 



(33) 
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Recently, in two investigations [0, |2^, effects of the transverse quark motion on 
diffractive charmonium production were discussed. In Ref. |^, the presence of a Q^ 



independent (!) correction was claimed, which contradicts the strict asymptotic QCD 



result of [0. For photoproduction, the correction term of Ref. |26] is by a factor of 24 
smaller than our leading twist, order 0{k^) correction of Eq. (^91) . To be able to compare 
with the result of Ref. |jl6|, we use the expression of T{Q'^) in transverse momentum 



space, i.e., Eq. (^. The correction factor for J/ip photoproduction discussed in Ref. 
T^ can be obtained from our T((5^ = 0) of Eq. (|^) by approximating At4>-y{z,kt) with 
the respective leading order expression in O {kf/m^) and by neglecting the longitudinal 
relative motion of the quarks, i.e., by setting (f)v{z^kt) = 6(z — |j (f)Y{kt). In addition, 
a Gaussian form for the wave function (f)v{kt) was assumed in Ref. [|T^]. All of these 
approximations diminish the relative contribution of large quark momenta, and hence 
result in a significantly weaker suppression. This was already pointed out in Ref. [Q in a 
footnote. 

2.6 The i-slope of diffractive vector meson production 

It was demonstrated in Ref. [0] that, in the limits of fixed small x and Q'^ —>■ oo, the 
t-slope of the vector meson electroproduction cross section should be flavor independent 
and determined solely by the slope of the gluon-nucleon scattering amplitude. However, 
the contribution of finite b ^ quark separations in the production amplitude of Eq. 
(|16|) - cf. Fig. 2 - lead to a Q^ and flavor dependence of the t-slope. This effect 
can easily be incorporated into Eq. ( [T6| ) by evaluating the matrix element of the factor 
g-izqt-b _|_ ^izqt-b _ ^-t{zqt+i)-b _ ^i{zqt+i)-b i^q^^qq^i the wavc fuuctious of the photou and the 

vector meson. Here, I is transverse momentum of one of the exchanged gluons, qt is 
the transverse component of the four-momentum transfered to the target nucleon and 
t = — l^tp (we neglect here terms proportional to x). This amplitude can be written, in 
factorized form, as a convolution integral over It of the hard blob and the non-diagonal 
gluon distribution in the target at a virtuality P. Since b is small, it is reasonable to 
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decompose this expression over powers of b up to b"^. After that, the respective amphtude 
factorizes into a product of the hard blob (accounting for the t dependence) and the non- 
diagonal gluon distribution in the target at virtuality ~ p-. Note that since the momentum 
integrals in the hard blob depend on 6^ only logarithmically, the b dependent term in the 
slope of the gluon-nucleon amplitude should decrease with b at least as 6^/ln^. So we 
can neglect it, to a first approximation, as compared to the effects of the form factor in 
the 'jI —^ V vertex, which are proportional to 6^. Besides, studies of soft elastic scattering 
indicate that even for such processes the main contribution to the t-dependence of the 
amplitude comes from hadron form factors (if the energy is not so large that Gribov 
diffusion contributes). Hence, we expect that, in the hard regime, the b dependent term 
in the amplitude will have a numerically small coefficient, in addition to being suppressed 
by the In -^ factor. 

Thus, effectively, one should include, similar as for a form factor of a QQ bound 
state [^, in the integral on the right hand side of Eq. ([T6| ) an additional factor of 
^-izqt-b^ where qt is the three-momentum transfered to the target nucleon and t = —\qt\'^- 
Parameterizing as usual ^ = Ae^^* for small t, we can calculate the Q^ dependence of 
ABviQ^) = BviQ^) - BviQ^ -^ oo), i.e., the contribution of the hard blob of Fig. 1 to 
the t-dependence of the cross section, from 

AB (O"^) = 1 Idz(Pbcf>^*^{z,b) (Tgg-jv(x,&) (^v{z,b) z^ b'^ 
2 Jdzd% (f)^*^{z,b) aqgN{x,b) (f)v{z,b) 

with the color dipole cross section crqqN{x,b) of Eq. (|TU|), and (f)^^{z,b) of Eq. (jT^) and 

4>v{z, b), the qq light-cone wave functions of the photon and the vector meson, respectively. 

Results of such a calculation are presented in Fig. 3 for J/\E' and p-meson production. 

Thus, the dependence of the cross section on t contains information on the distribution of 

color in the produced vector mesons. Note that the experimentally observed t-slope for 

J/\E' photo- and electroproduction and p" electroproduction at large Q'^ is of the order of 

By ~ 4 — 5 GeV^^. The main conclusion from Fig. 3 is thus that the t-slope of diffractive 

vector meson production is determined mostly by the gluon-nucleon scattering amplitude, 

the differences in By between different flavors are small for realistic Q^, and they vanish 
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in the Q^ ^ oo limit. 
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Figure 3: The contribution of the hard blob to the t-slope of diffractive electroproduction 
of longitudinally polarized p° and J/ip vector mesons. 

We have demonstrated in Ref. 0] that, at sufficiently small x, i.e., close enough to the 
low a;-range probed at HERA, higher twist effects may become important. This would 
also lead to a break-down of the universality of the t-slope. This effect can be estimated 
by including double scatterings of the qq pair off the nucleon |3^, H, |3ll . Neglecting the 
small difference of the average b for single and double scattering, we can calculate the 
i-slope of the rescattering amplitude from 



da 
~dt 



da 
It 



t=o 



.Bt/2 



1 {(^qqA^.b)) 



(e^*/4 + re^'/' 



Here ^ 



t=o 



16nB {aggNix,b)) 
is the cross section given in Eq. (j^) and 

(a2-^(x,6)) _ Jdzd'^b(f)^i{z,b) a^g^ix^b) (f)v{z,b) 



(35) 



(36) 



{aqqNix.b)) Jdzd% (j)Y^{z,b) aqqNix,b) (j)viz,b) 

where aqq^^x, b) is again the color dipole cross section of Eq. (|r^). The quantity denoted 
r is the ratio of the inelastic to elastic diffractive production of vector mesons at large 
Q^, and experimentally r ~ 0.2. B is the slope of the inelastic production, i.e., the 
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7l + P ~^ V + X process, and it is, so far, not well known experimentally. Since for 
this process, in difference from the elastic production, there is essentially no form factor 
at the nucleon vertex, the quantity B is much smaller than the elastic slope B, and a 
natural guess is S ~ 1 — 2 GeV~^. Experimentally [IS], the ratio of inelastic to elastic 
J/\I' production is of the order of 0.5 — 0.7, i.e., ^ ~ 0.5 — 0.7. 

For our numerical estimates we set r = 0.25, B{x ~ 10^^) = 5 GeV^^, and B = 2 
GeV^^. In Fig. 4a we show the t-dependence of the diffractive p° electroproduction cross 



section, i.e.. 



da 
dt 



I 



da 
dt 



t=0 



of Eq. 



Since the color dipole cross section <JqqN{x,b) 



of Eq. ([T0|) is proportional to xGn{x, X/b'^) oc a;0-2-o.3^ gq_ ( psf ) leads to an increase of 
the t-slope with decreasing x. This can be seen from Fig. 4a, where we compare ^ 
for X = 10"^ (dashed line) and x = 10~^ (solid line) with the leading twist result e~^* 
(dotted line). 
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Figure 4: The t-dependence of the diffractive vector meson production cross section, i.e.. 



da 
dt 



/ 



da 
dt 



t=0 



of Eq. (|35[), and the change of the t-slope with energy, i.e., a'(t) of Eq. 



(P7|). Results are shown for p" electroproduction at Q^ = 10 GeV^. 

The change of the t-slope with "energy" s = — is usually parameterized in the form 



B{s) = B{so) + 2a'{t) In 



s 
So 



(37) 



and the quantity a' increases with —t. This can be seen from Fig. 4b, where we show a'{t) 
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as a function of t for various x. Again, due to the increase of the color dipole cross section 
CTqqi\i{x,b) with energy, the increase of a' with —t is more dramatic for smaller x (larger 
energies). Note that the numerical results shown in Fig. 4 refer to p° electroproduction 
at Q^ = 10 GeV^. As can be seen from Fig. 2, the respective Qlft is very similar to the 
Qlff relevant for J/ip photoproduction, and hence the numerical estimates shown in Fig. 
4 should thus be approximately valid also for J/ip photoproduction. Fig. 4 suggests that 
a study of the t-slopes of diffractive vector meson production may yet provide another 
sensitive probe of the dynamics of hard diffraction. 

3 The quarkonium light- cone wave function 

In order to be able to evaluate the asymptotic correction rjv of Eq. (H) as well as the 
T(Q2) and R{Q^) of Eqs. (0) and (g) or (|2|) and (^, we need the light-cone wave 



function of the qq leading Fock state in the vector meson. We will discuss this quantity 
in detail in this section. Note also that, as a result of the factorization theorem in QCD, 
it is the distribution of bare (current) quarks that enter in the description of these hard 
processes, and therefore, a priori, there should be no simple relation between this quantity 
and non-relativistic potential models. 

3.1 Non-relativistic potential models 

Due to the large value of the quark mass, it is generally assumed that a non-relativistic 
ansatz with a Schroedinger equation and an effective confining potential yields a fairly 
good description of heavy quarkonium bound states. The various models - see Ref. ^^ for 
an overview - differ in the functional form of the potential, but they all give a reasonable 
account of the cc and bb bound state spectra and decay widths. The same holds for the 
light-front QCD bound state calculation of Ref. [@. In Fig. 5, we display the quantities 
Roo{r) (normalized such that Jdrr"^ |i?oo('")P = 1) and 47r^A;^0v(fc) (normalized such that 
■7^^ Jd^k \(f)v{k)\ = 1). For the latter, we also plot a Gaussian fit adjusted to reproduce 

21 



0y(fc) at small k. It turns out that the wave functions can be well approximated at small 
k by Gaussians, while, at large k, they fall off much slower and they display a significant 
high momentum tail. 
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Figure 5: The non-relativistic quarkonium wave functions for the heavy ground state 
mesons J/ip and T from various potential models [|, ||, 0, |[ and a light-front QCD bound 
state calculation [^. In the lower part of the figure, we also show a Gaussian fit adjusted 
to reproduce 0y(/c) at small k (dotted lines). 

Note that for our actual numerical calculations we will restrict our considerations to 
the models of Refs. 0, and §] for which the mass of the constituent quark is close 
to the mass of the bare current quark, i.e., rric ~ 1.5 GeV and mf, ^ 5.0 GeV. This is 
necessary to keep a minimal correspondence with the QCD formulae for hard processes 
which are expressed through the distribution of bare quarks Q. 
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In Fig. 6, we show the contributions of the different regions in momentum space to 
the integral Jd^k(j)v{k) for the potential model of Ref. P] (logarithmic potential). This 
integral appears, for instance, in the expression for the V -^ e~^e~ decay width. Especially 
for J/ip mesons, the conventional non-relativistic potential models lead to a significant high 
momentum tail in the respective wave functions, and the contribution of the relativistic 
region^ - > 1 (or k > m) to the integral Jd^k(f)y{k) (the shaded area in Fig. 6) becomes 
large. For the potential model of Ref. |^, the contribution of the relativistic region k > m 
to the integral under consideration is 30% for J ftp (and < 10% for T). Also, for the 
J/ip, half of the integral comes from the region k > 0.7m. This is in line with the QCD 
prediction of large relativistic corrections to the cc bound state equations [^ , and it puts 
the feasibility of a non-relativistic description of heavy quarkonium production in high 
energy processes seriously into question. 
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Figure 6: Histogram of the relative contributions of the different regions in momentum 
space to the integral J(fk(j)v{k) for the potential model of Ref. p. 

The fact that, in particular for the J/ip meson, our numerical analysis yields a sig- 
nificant value for the high momentum component in the respective non-relativistic wave 
functions is a very important result which should have consequences far beyond the topic 
of diffractive vector meson production. The large high momentum tail, visible in the lower 
^Evidently, relativistic effects should become important already at significantly smaller k. 
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part of Fig. 5, and the significant contribution of the relativistic region to the integral 
Jd^k 0y (fc), displayed in Fig. 6, indicate that the J/ip meson is not really a non-relativistic 
system! This puts the non-relativistic ansatz employed in the various potential models 
@! i) 0! § as well as in the light-front QCD bound state calculation seriously into 
question. 

However, there are more inconsistencies between the non-relativistic ansatz and the 
hard reaction considered here. For once, the requirement of self-consistency dictates that 
since in our formulae we use the gluon distribution xG'iv(x,(5^) extracted from the data 



within a certain renormalization scheme (MS), we are indebted to use the bare quark 
mass defined within the same scheme. This means that, in our final formulae in Eqs. (|^), 
(0)' (H)' (P^; (P^, (P7|) and (pSD, the pole or constituent quark mass m has to be replaced 
by the running mass rriruniQlff), where ||3^ 

3a, 
3^ 



m' - mUQlff) = m^(l-^) . (3^ 



Here, as is evaluated at Qljf, i.e., at the effective scale of the reaction determined via 
the so-called "rescaling of hard processes". This is another consequence of the differ- 
ence between soft, non-perturbative physics (as described, for instance, by non-relativistic 
quarkonium potential models) and hard perturbative QCD, and it further stresses the 
inadequacies of a naive straightforward application of non-relativistic potential models in 
this context. 

Another mismatch between the soft non-relativistic and the hard light-cone approach 
appears within the evaluation of the V — > e+e^ decay width. When Fy^e+e- is calculated 
from the non-relativistic wave function (f)v{k) via 



V^e+e- ,.2 



levra^e^ /^ 16a, 



M^ V Sir 



Mk) 



(39) 



(27r 

a QCD correction factor |12| appears, 1 — ^|^, which can be numerically large (^|f^ ~ 
0.35—0.65 for J/ip), while no such term is present within the relation p[ with the light-cone 
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qq wave function 4>v{z, k, 



t), 



V^e+e- 



Mv 



dz 



(fkt 
167r3 



<Pv{z,kt) 



(40) 



The appearance of this correction factor is the main differences between the various non- 
relativistic potential models and a "true" QCD approach, in which the light-cone wave 
function of the minimal qq Fock component in the vector meson is employed. It is a radia- 
tive correction to the matrix element of the electromagnetic current calculated, essentially, 
while neglecting quark Fermi motion effects. The respective Feynman diagram is shown 
in Fig. 7. 
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Figure 7: The QCD radiative correction [12| to the V — >■ e e decay width 



The correction arises from the exchange of a gluon between the quark and the antiquark 
in the vector meson with fairly large transverse momentum, (It) ~ m. The physical 
interpretation of the 1 — i|^ correction factor is that it "undresses" the constituent quarks, 
which are the relevant degrees of freedom of the non-relativistic wave function, back to 
current quarks, which, in turn, are the degrees of freedom the light-cone wave function 
refers to and to which the V -^ e"*"e~ decay width is ultimately connected. This, once 
more, underlines the limits of applicability of non-relativistic potential models in that 
context. Note that this radiative correction is also present in the light-front QCD bound 
state calculation of Ref. |^ because also there the relevant degrees of freedom are dressed 
constituent and not bare current quarks. 
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3.2 The light-cone wave function 

Leaving these issues behind for the moment, we can, in principle, deduce a hght-cone 
wave function 0y(-2, kt) appropriate for the evaluation of time-ordered perturbation theory 
diagrams from the non-relativistic wave function (f)v{k). This requires a translation of 
conventional non-relativistic diagrams into light-cone perturbation theory diagrams. This, 
in turn, can be achieved by the purely kinematical identification of the Sudakov variable 
z, which denotes the fraction of the plus-component of the meson's momentum carried by 
one of the quarks, with 

- = ifi + ^fi=?V (41) 



2 V Vt 



2 _|_ ^^2 



This yields 



k^ -. ''^ + ('" - '^'"^^ (42) 

" 4z{l-z) ' ^ ^ 



Jkf + 711? 
d'k -^ i[^(rr^jp^-^^*' (43) 

where ±kt are the quarks' transverse momenta. This, together with the conservation of 
the overall normalization of the wave function, 

then gives a relationship between the light-cone and the non-relativistic wave function: 
(f)viz,kt) ■ 



kl + m? 

V k 



kl + {2z - Ifm^ , ^^^^ 



\A[z{l-z)Y \ \ Az{l-z) 

From (f)v{z,kt) we then calculate the quarkonium's wave function in transverse impact 
parameter space (j)v{z,h) through a two-dimensional Fourier transformation, 

MzM = l^e^'^-'Mz.h). (46) 

Obviously, the non-relativistic quarkonium model, designed as a description of the qq 
constituent quark component - including the 1 — i|^ factor which accounts for radiative 
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corrections - does not include gluon emission at a higher resolution. So it is not surprising 
that the (pv{z,b) that we find does not display the expected asymptotic behavior ||TT| , 



Mz,b = 0) 



zil 



(47) 



This is illustrated in Fig. 8. There, we compare the quarkonium wave functions ipviz-, b = 
0) obtained in that manner from the non-relativistic potential models of Refs. 0, 
and [§] with a hard wave function 0y"'''^(z, 6 = 0) = aoz{l — z), where the parameter oq 
was adjusted by means of Eq. ( ^OD to reproduce the vector meson's leptonic decay width 
ry_>e+e-. For trausversely polarized vector mesons, the light-cone wave function should 
behave as oc z^(l — 2;)^. This follows from the analysis of respective pQCD diagrams. 
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Figure 8: The quarkonium wave functions, 0j/^(z, 6) and <^-x{z,\)), for 6 = 0. The dot- 
dashed, dashed and dotted lines correspond to the non-relativistic potential models of 
Refs. 0, p and 0, respectively, and the solid fines refer to tfie "fiard pfiysics" limit 

0y(2;, 6 = 0) oc z{\ — 2;). 

3.3 Hard physics 

We argued in subsection ^?T] tfiat, altfiougfi tfie average properties of fieavy quarkonium 
bound states migfit, in general, be quite well described witfiin a non-relativistic framework 
due to tfie large value of tfie quark mass, tfiere appear significant fiigfi momentum (or 
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"hard physics" ) corrections if observables are considered which crucially depend on short 
distances. 

One example is the leptonic decay width, Ty^e+e-, which acquires large radiative 
corrections in a non-relativistic potential model. An analysis of the respective Feynman 
diagram, shown in Fig. 7, yields that these corrections arise from relativistic momenta, k ^ 
m. Even putting those corrections aside, already the quantity which is related to the decay 
width in zeroth order, Jd^k(f)v{k), contains large contributions from the relativistic region 
k > m (30% for J/ip for a typical potential model). And when we (purely kinematicly) 
translate the non-relativistic wave functions into light-cone coordinates, we find that they 
do not display the expected asymptotic short distance behavior (f)y{z,b = 0) oc z{l — z) 
as dictated by perturbative one-gluon exchange. 

This suggests that the non-relativistic potential model wave functions might describe 
the qq leading Fock state in heavy quarkonia for fairly large (average) distances, but the 
description breaks down in the limit of small distances or high momenta. As these play a 
crucial role for the processes we are interested in, we designed the following strategy: 

First, we extract a light-cone wave function from a non-relativistic potential model 
through the purely kinematical transformations of Eqs. (^TJ) through (^5]), which we then 
Fourier transform into transverse impact parameter space via Eq. (^61). However, we have 
confidence in that wave function, which we denote (j)y^{z, b), only for transverse distances 
6 > — , and we expect it to be modified at shorter distances by means of the "hard physics" 
corrections discussed in the above. We thus set 

r 0^^(2,6) for6>6o, 
Mz,b) = r^ ; '; - ' (48) 

where fen ~ — • 

" m 

The wave function (py"{z, b) is then constructed such that: 1) (pv{z, b) and d(t)v{z, b)/db 
are continuous at 6 = b^, 2) 4>y^ {z, b) has the correct asymptotic behavior dictated by the 
perturbative exchange of hard gluons, i.e., (py^^z^b = 0) oc z{l — z), and 3) (py^i^z.b) 
reproduces the vector meson's leptonic decay width without account of the radiative cor- 
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rection term 1 — ^^^, i.e., Eq. ( ^OD is used to calculate Ty^e+e-- We expand 4>v"{z, h) in 
terms of Gegenbauer polynomials^, 

0t^^(z,6) = a,{h)z{l-z) (l + E an{h)Cl'\2z-l)\ , (49) 

V n=2,4,- / 

and we assume that the coefficients aj(6) depend on 6^ only through second order, i.e., 
aiip) = ttjo + anb"^ + ai2b'^. This, together with the conditions 1) through 3) is sufficient to 
unambiguously determine (f)y^{z,b). In our actual numerical calculations, we set bo = 0.3 
fm for J/ip and 6o = 0.1 fm for T. Respective wave functions are shown in Fig. 9. The dot- 
dashed, dashed and dotted lines show the non-relativistic wave functions (py^{z = |,6) 
before the "hard physics" corrections discussed in this subsection were imposed, and the 
solid lines depict the modified wave functions (py'^Xz = |, 6) of Eq. (|49D. 

Note that the "hard physics" corrections, which we introduced in the above, address 
effects that are of higher order in an expansion in — . But the prescription of modifying 
the wave function at 6 < 6o only accounts for some (but not all) corrections to this order. 
We thus emphasize that the corrections outlined in this subsection are a model estimate 
only. 

3.4 Vector meson production 

For the potential models of Refs. 0, [Q and |^, the non-relativistic wave functions 
(f)Y^{z, b) yield values for the asymptotic correction factor rjv of Eq. (|^) of rjj/^ ^ 2.3 — 2.4 
and 7]x ~ 2.1 — 2.2 if the "hard physics" correction outlined in the last section is not 
considered. While, with that correction, they yield r]v = 3. Note that the static limit, 
i.e., (pviz^b) = 5{z - l/2)0y(6), gives r]v = 2. 



Furthermore, in line with the discussion in subsection 3.1, we do not use the pole 



mass m in our final formulas, but we replace it with the running mass rrirun, as given 
■^The expansion in Gegenbauer polynomials has nothing to do with renormalization group methods. 
They provide a complete basis for the (J3viz,b) under consideration and allow a smooth interpolation 
between the & — > and b ^ l/m regimes. The series in Eq. (H9) is terminated when convergence is 
achieved, which, in practice, is at n ^ 10. 
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Figure 9: The quarkonium wave functions, (f)j/^{z,b) and (j)r{z,b), for -2 = |. The dot- 
dashed, dashed and dotted hues correspond to the non-relativistic potential models of 
Refs. 0, and 0, and the solid lines refer to the inclusion of the "hard physics" 



corrections of Eqs. (^H]) and (^) for b < bo. We set bo = 0.3 fm for J/ip and bo = 0.1 fm 
for T. 



by Eq. (^Sp. We can then use the wave functions 0y(z,6), that we constructed in the 
last two subsections, to calculate the correction factors of Eqs. (^ and (H) or (p5|) and 
(p6|). Putting everything together, we can rewrite the forward differential cross section 
for photo- and electroproduction of heavy vector mesons of Eq. (^ as the product of an 
asymptotic expression and a finite Q^ correction, C((5^), where 



dt 



t=o 






(^siQlff) (1 + iP) xGn{x, Qlff] 



1 + e 






with 



CiQ' 



wY (_Q 



2 I 4^2 \4 



3 J \Q^ + Amj, 






1 + e 



C{Q') ,(50) 



(51) 



m 



Here, rjv is the leading twist correction of Eq. (Q), the factor T^Q"^) of Eq. (0) accounts for 
effects related to the quark motion in the produced vector meson, e is the (virtual) photon's 
polarization, and the factor R{Q'^) of Eq. (Rl) parameterizes the relative contribution of 
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the transverse polarization. The pole mass m we set to m = 1.5 GeV for J/ip and to 
m = 5.0 GeV for T production, and rrirun is the "running mass" of Eq. (^) which, 
through Qlff, depends on Q^ and the vector meson's wave function. 

Results for the Fermi motion suppression factor, T{Q'^) of Eq. (^, and the finite 



Q^ correction, C(Q^) of Eq. (0), are shown in Fig. 10 for J/ip and T photo- and 
electroproduction. The calculations are based on vector meson wave functions from the 
models of Refs. ||^ and 0. The solid line, labeled hard, refers to the inclusion of 
the "hard physics" corrections of Sect. |3.3| . For the evaluation of C((5^), the photon's 
polarization e was set to 1. 
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Figure 10: The Fermi motion suppression factor, T(Q^) of Eq. 
correction, C{Q'^) of Eq. (|5T|), for J/ip and T production. 



and the finite Q^ 
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It can be seen from that figure that, for reasonable Q^, the correction factor C(Q^), 
which measures the suppression of the cross section due to the quark motion in the pro- 
duced vector meson, is significantly smaller than 1. This shows that the asymptotic ex- 
pression, i.e., Eq. (0) with the finite Q^ correction C((5^) set to 1, is valid for extremely 



large Q^ only. Note that the "hard physics" corrections of Sect. |3.3| lead to a stronger 
suppression in T((5^), but, at least for J/-?/' production, to less suppression in the final 
correction factor C[Q'^\ The reason for this is, firstly, that the "hard physics" correction 
increases r/y of Eq. (^ from around 2.1 — 2.4 to 3, and, secondly, that the running mass 
mrun of Eq- (|38D is smaller than the pole mass, which also enhances the cross section. In 
addition, the relative contribution of the transverse polarizations R{Q^) of Eq. (^) is very 
close to 1 both for J/ip and T production for all experimentally accessible Q^ if the "hard 
physics" corrections are left out. However, at least for J/ip production, after the "hard 
physics" corrections are considered, R{Q^) increases significantly with Q^. This, together 
with the changes through rjv and rrirun lead to the difference between T{Q'^) and C{Q^). 
The cross sections are enhanced also due to the so-called "rescaling of hard processes" , 
because the virtuality that enters in the gluon density, Qlff of Eq. (p^, is larger than 
the naive estimate oi Q = — ■^—^. This was discussed in detail in Section [27^ . 

Note that for photoproduction, i.e., for Q^=0, only the transverse polarizations are 
present, and the correction C^Q"^) of Eq. ( |5TD takes on the form 



C(0) oc T(0)i?(0) oc 



I z^t,)2 Jd'^h (j)v{z,kt) At (j).,{z,kt] 

dz 



;nhances smearing in th 
i.e., the contribution of asymmetric qq pairs with z ^ \ is pronounced. 



(52) 



The presence of the ^2^-^-12 fei"ni strongly enhances smearing in the longitudinal motion. 



3.5 The ratio of T and J/ijj photoproduction 

One can furthermore conclude from Fig. 10, together with our master formula in Eqs. 
(0) and (|5lD , that, after an eventual luminosity upgrade, a significant production of T 
mesons is expected at HERA. The cross section ratio of T to J/ip photoproduction (at 
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fixed x) is approximately 

a{^ + p^ T +p) _ rxMlm^ |a,(l + z/3)xG^(Q2^/T])|2 Ct(0) 



• (53) 



a(7 + p^W+p) Tj/^Mj/^ml |a,(l + z/?)xG^(Q2^/J/^])|2 Cj/v,(0) 

The first factor on tlie riglit liand side of Eq. (^) is tlie dimensional estimate, and it 
yields a relative suppression of T photoproduction of about 1 : 2000 if we set for the 
quark masses rric = 1.5 GeV and rrib = 5.0 GeV. The second term arises due to the so- 
called "rescaling of hard processes", and it enhances the cross section ratio by a factor 
of about 3 for x = 10"'^. The third term is connected to the wave function dependent 
effects, and it enhances the production ratio also by a factor of about 3. All together, the 
cross section for T photoproduction is suppressed by approximately 1 : 200 as compared 
to J/ip photoproduction for the same x. For the same W, an extra suppression factor 
^ (Mr/Mj/^)^'^ ^ 2.4 is present. Note that the different Q^ scale and higher twist effects 
(the "rescaling of hard processes" as well as the C((5^) correction) increase the relative 
yield that we predict by about an order of magnitude as compared to the naive dimensional 
estimate! 

4 The J/^ photo- and electroproduction cross section 

In Fig. 11, we compare our prediction^ for the J/\E' photoproduction cross section with 
the data. We used a slope parameter of Bj/,^ = 3.8 GeV^^, as measured by the HI 



collaboration |T^ , to calculate the total cross section from our predictions for the forward 
differential cross section at t = 0, and the Fermi motion corrections and the "rescaling 
of hard processes" are accounted for. For the former, the charmonium potential of Ref. 
was employed and the "hard physics" corrections, as outhned in Sect. |3^ , were taken 
into account. We furthermore replaced the quark pole mass with the running mass rrirun 
from Eq. (|38|), and we set x = — W2-^- The formulas to obtain the forward differential 
cross section are given in Eqs. (|50| ) and (pl). 

''A respective FORTRAN program is available by request from koepf@mps.ohio-state.edu or via the 



WWW at http://www.physics.ohio-state.edu/~koepi. 
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Figure 11: The J/\l/ photoproduction cross section for several recent parameterizations of 
the gluon density [^ |35| , |36| in comparison with experimental data from E401 [0, E516 
[, E687 pg, ZEUS '93 [M, and HI |T7|. 



As can be seen from Fig. 11, the predictions of our pQCD calculation agree with the 
data within the uncertainties in the nucleon's gluon density, and the energy dependence 
of the data is much better reproduced within the pQCD picture, where a oc p^o-^-o-s^ 
than through the soft Pomeron model [^, where o oc W^''^'^ . A rough fit |T^ to the data 
depicted in Fig. 11 yields a oc iy°'^. 

To investigate the Q^ dependence of J/\E' production, we show in Fig. 12 the ratio of 



the electro- to photoproduction cross sections, i.e., we plot — 



»+p->.j/i/»+p(Q ) 



^/t2— ^ as a function 

0"7+p->J/i/'+pW — uj 

of the virtuality of the photon. In particular, we compare a calculation where the Fermi 
motion corrections were left out (dotted lines), with an evaluation were the latter effects 
were included while either using just the non-relativistic wave functions (dashed lines) or 
also accounting for the "hard physics" corrections (solid lines). 

As can be seen from Fig. 12, the Fermi corrections are necessary to achieve agreement 
with the data. However, at this point, the quality of the data is not sufficient to distinguish 
between the various potential models or to decide whether the "hard physics" corrections 
which were imposed on those wave functions at small transverse inter-quark distances - 
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Figure 12: The ratio of the J/\l/ electro- to photoproduction cross section for two recent 
parameterizations of the gluon density p4| , |35| and for various potential models P, ^, ^ 
in comparison with experimental data from HI |r?|. 



see Sect. p.3| - lead to an improvement. This should change if the '95 data, which have 
much better statistics, become available. The fact that we somewhat underestimate the 
J/ip photoproduction cross section - see Fig. 11 - and, at the same time, overestimate the 
Q^ dependence of J/ip electroproduction - see Fig. 12 - suggests that our quark motion 
correction factor C{Q'^) of Eq. (|5T1) is too small at Q^ = and it falls off too quickly 
at larger Q^ . This implies that the wave functions which we use fall off too slowly in 
transverse momentum space and they are too steep as a function of the impact parameter 
6, i.e., the respective (fc^) is too large. 



5 The p^ electroproduction cross section 

Although the main topic of this work is heavy meson photo- and electroproduction, we 
still consider an update of our predictions of Ref. 0] in regards to p° electroproduction 
warranted in light of the new data as well as theoretical developments in that realm. 
Currently, absolute cross sections for exclusive p-meson production are available from 
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NMC |^T|, ZEUS 1^, and HI |T8[, and preliminary results exist from ZEUS from the 



1994 run [^ . From our predictions for the forward differential cross section, , 

i=0 

of Eq. (Q), the total cross section was calculated using a slope parameter oi Bp = 5 GeV~^. 
This is consistent with the values given by the NMC ^ (4.6 ± 0.8 GeV"^) and ZEUS 



(5.1 ±1.2 GeV ^) collaborations and slightly smaller than that obtained by HI |T8 
(7.0 ±0.8 GeV-2). 

To indicate separately the spread that arises from the different available gluon densities 
and the uncertainty that stems from the various proposed p-meson wave functions, the 
theoretical predictions are shown for two (extremal) gluon densities, GRV94(HO) of Ref. 



3^ and MRSR2 of Ref. |^5|], and two different wave functions, termed "soft" and "hard". 
The "soft" wave function refers to a (f)p{z, kf) oc expf — ^.-^ '^. j with an average transverse 
quark momentum of {k'f) = 0.18 GeV^ as extracted from a QCD sum rule analysis by 
Halperin and Zhitnitsky |2^, and the "hard" wave function corresponds to a (j)p{z,kt) oc 



z{l — z) ,^'1 2p obtained in another QCD sum rule analysis (for pions) by Lee, Hatsuda 
and Miller ^^. For the latter, (fc^) = 0.09 GeV^. As outlined in detail in the above, 
the wave function enters through the Fermi motion suppression factor, T(Q^) of Eq. (^). 
T{Q^) is depicted in Fig. 13 for various available p-meson wave functions: "hard" and 
"soft" were discussed in the above, "softi" refers to a duality wave function of the form 
Of So — ^(i'Lz) ) with (kf) = 0.15 GeV^ obtained in Ref. |4^ and used for a similar analysis 
in [^, and "soft2" labels a two-peak Gaussian favored in the analysis of Halperin and 



Zhitnitsky pO|. Note that the latter wave function seems quite extreme as it would 



correspond to a transverse spread of the qq component which is larger than the meson's 
size! 

The comparison of our prediction^ with the most recent experimental data is shown 
in Fig. 14. In the kinematic domain were our approach is expected to be applicable, 
X ^ 0.01 and/or W ^ 30 GeV, our predictions agree with the data within the spread 



^A respective FORTRAN program is available by request from koepf@mps.ohio-state.edu or via the 



WWW at http://www.physics.ohio-state.edu/~koepi. 
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Figure 13: The Fermi motion suppression factor T{Q^) of Eq. (^ for p° electroproduction 
for various p-meson wave functions from Refs. |^, [Q and ||4^ . 



through the various parameterizations for the gluon density and the uncertainty which 
stems from the vector meson's wave function. Note, in particular, that as Q^ increases 
the energy dependence of the prehminary ZEUS data |T^ approaches more and more the 
hard physics hmit, a oc ^yO'^'-O-S^ which is very different form the soft Pomeron prediction 
||l|], o" oc p^o-22-o.32_ ^jj-^jg could indicate a transition form soft to hard physics in the Q^ 
range depicted in Fig. 14. 

There are two reasons why our predictions should not really reproduce the data very 
well at smaller Q^. Firstly, smaller Q^ correspond to larger transverse distances, and hence 
the pQCD approach outlined here loses some of its validity. Secondly, at very small x, the 
increase of these cross sections with energy is restricted by the unitarity of the S'-matrix, 
and even more stringent restrictions follow from the condition that the leading twist term 
should be significantly larger than the next to leading twist term [Q. The kinematical 
region where this limit becomes important moves to larger x for decreasing Q^. However, 
whether the softer energy dependence of the cross sections at the smaller Q^ is really due 
to the unitarity limit slow-down is unclear at the moment. Further work is in progress in 
that realm E5 . 
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wave functions in comparison with preliminary ZEUS data [19]. 
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The Q^ dependence of the cross section is commonly parameterized through a quantity 
a, where (for fixed W^ 

a{-f* + p -* p" + p) oc g-2" , . (54) 



The various experiments yield a = 2.1 ± 0.4 |42|, a = 2.4 ± 0.3 [19| and a = 2.5 ± 0.5 



T8| at {Q"^) ^ 12 GeV^ and (W) ~ 80 GeV. Neglecting the Fermi motion corrections 



and the "rescaling of hard processes", our theoretical predictions yield a ~ 3.3 without 
the corrections, while we find a ~ 2.6 if we take the quark motion and "rescaling of 
hard processes" into account. To evaluate the correction factor T{Q'^) of Eq. (p5D , we 
again used the wave function (l)p{z,kt) oc expf — ^^_'^. j with an average transverse quark 
momentum of (kf) = 0.18 GeV^ as extracted from a QCD sum rule analysis |2^. Hence, 



our predictions agree with the measurements only if the Fermi motion corrections and 
the "rescaling of hard processes" are taken into account. This underlines our claim |Q 
that the Q^ dependence of those cross sections could eventually be used to probe the 
transverse momentum distributions within the produced vector mesons. However, at 
present, the data are still far too crude to extract conclusive information on this quantity. 
Note, furthermore, that our prediction refers to the Q^ dependence of the longitudinal 
cross section, ai, while the experimental values listed in the above correspond to the Q^ 
dependence of the total cross sections, a = ar + eai. 

6 Conclusions 

In this work, we focused the QCD analysis of Refs. and [Q on heavy quarkonium {J/ip 
and T) photo- and electroproduction, and we extended the respective formalism, which in 
Refs. PI and 0] was applied to the production of longitudinally polarized vector mesons 
only, to transverse polarizations as well. 

For non-asymptotic momentum transfers, the respective hard amplitude is sensitive 
to the transverse momentum distribution in the qq light-cone wave function of the lead- 
ing Fock component in the produced vector meson. This leads to a suppression of the 
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asymptotic predictions, i.e., to an interplay between the quark(antiquark) momentum dis- 
tribution in the vector meson and the Q^ dependence of the corresponding cross section. 
We derived the respective expressions for the Fermi motion suppression factor, T{Q'^) of 
Eqs. (^ and (pSj), and the relative enhancement of the transverse cross section, R{Q^) of 
Eqs. (D and (|2|), to leading order in ^2+4^2 • 

The evaluation of these factors required a detailed study of the vector meson's qq light- 
cone wave function. Motivated by the large value of the quark mass in heavy quarkonia, we 
started from conventional non-relativistic potential models, which we critically examined 
and confronted with QCD expectations. In particular for the J/ip meson, our numerical 
analysis yields a significant value for the high momentum component in the respective 
wave functions, visible in the lower part of Fig. 5, and a significant contribution of the 
"relativistic region" - > 1 to the integral Jd^k(j)v{k), displayed in Fig. 6. This is in line 
with large relativistic corrections to the corresponding bound state equations [|1^. These 
large relativistic effects question the feasibility of a description of heavy meson production 
in high energy processes based on a non-relativistic ansatz. This is a very important 
result which should have consequences far beyond the scope of diffractive vector meson 
production, and it indicates that the J/ip meson is not really a non-relativistic system! 
We therefore designed an interpolation for the wave function of heavy quarkonia which 
smoothly matches the results obtained from non-relativistic potential models with QCD 
predictions at short distances. 

We then used the latter to evaluate the finite Q^ corrections for diffractive J/ip as well 
as T production. We find fairly good agreement of our predictions with the J ftp data, and 
we predict a measurable production of T mesons at HERA - especially after a luminosity 
upgrade. We also update our comparison of longitudinal p° electroproduction with the 
data, putting special emphasis on preliminary ZEUS '94 data [jl9| that became available 
only recently. 

The discussion in this work affirms that hard diffractive vector meson production is 
exactly calculable in QCD in the same sense as leading twist deep inelastic processes. This 
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holds if only short distances contribute, which is the case for heavy flavors or production of 
longitudinally polarized p° at large Q^. The respective amplitude is expressed through the 
distribution of bare quarks in the vector meson and the gluon distribution in the target. 
This is qualitatively different from an application of the constituent quark model to these 
processes, as in Refs. [|l5l and [|l^. On the other hand, it makes these processes an ideal 



laboratory to study the qq leading Fock state in vector mesons. 
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